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Collective flavor oscillations are known to bring multiple splits in the supernova (SN) neutrino and 
antineutrino spectra. These spectral splits depend not only on the mass hierarchy of the neutrinos 
but also on the initial relative flux composition. Observation of spectral splits in a future galactic 
supernova signal is expected to throw light on the mass hierarchy pattern of the neutrinos. However, 
since the Diffuse Supernova Neutrino Background (DSNB) comprises of a superposition of neutrino 
fluxes from all past supernovae, and since different supernovae are expected to have slightly different 
initial fluxes, it is pertinent to check if the hierarchy dependent signature of collective oscillations 
can survive this averaging of the flux spectra. Since the actual distribution of SN with initial relative 
flux spectra of the neutrinos and antineutrinos is unknown, we assume a log-normal distribution for 
them. We study the dependence of the hierarchy sensitivity to the mean and variance of the log- 
normal distribution function. We find that the hierarchy sensitivity depends crucially on the mean 
value of the relative initial luminosity. The effect of the width is to reduce the hierarchy sensitivity 
for all values of the mean initial relative luminosity. We find that in the very small mixing angle 
(613) limit considering only statistical errors even for very moderate values of variance, there is 
almost no detectable hierarchy sensitivity if the mean relative luminosities of u^. and Ue are greater 
than 1. 
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I. INTRODUCTION 

Neutrinos coming from supernova explosions can give 
rich information about the explosion mechanism as they 
are the only particles that come from regions deep inside 
the core. They are also extremely useful for determining 
neutrino properties as neutrinos from SN1987A have am- 
ply demonstrated jJ-Q- In particular, it has been shown 
that the neutrino mass hierarchy and smallness of ^13 
can be probed using SN neutrinos 0, Q • In the last few 
years, focus has been on the effect of neutrino-neutrino 
interaction in the central regions of the core of supernova, 
giving rise to the so-called collective effects |6l-|43| . 

The 'collective' nature of simultaneous flavor conver- 
sions of both neutrinos and antineutrinos give rise to 
'spUts' in the spectra of the neutrinos and antineutrinos. 
These splits occur due to sudden change in the oscillation 
probability, causing spectral swaps which may end up in 
observable effects. Interestingly the impact of collective 
oscillations on the spectra are different for the Normal 
Hierarchy (NH) and the Inverted Hierarchy (IH). This 
opens up the possibility of identifying the neutrino mass 
hierarchy via observation of collective effects in the neu- 
trino signal from a future galactic supernova event (soj . 

However with the very small rate of occurrence of 
galactic supernova events (a few per century) one is 
forced to think of strategies of detecting the above men- 
tioned effect otherwise. One of the promising possibili- 
ties is the detection of Diffuse Supernova Neutrino Back- 
ground (DSNB) in the near future [ii - fiet . The cumu- 
lative number of neutrinos and antineutrinos produced 
by all earlier SN events in the universe result in a cosmic 
background known as DSNB. Though these neutrinos are 
yet to be detected, one has reasons to believe that it may 
be possible to observe them in the near future [47h66| . 



Present day upper limits are 1.2 cm~^ s~^ with ener- 
gies above 19.3 MeV for Super-Kamiokande (SK) water 
Cerenkov detectors [g^) and 6.8 xlO^ cm~^ s~^ with 
energies between 25 MeV and 50 MeV for Liquid Scintil- 
lator Detector (LSD) at 90% C.L. for both. 

The two key ingredients in the calculation of DSNB are 
(i) the SN rate which is proportional to cosmic star for- 
mation rate and (ii) the v and V energy spectra. Whereas 
reliable estimates are now available for the star forma- 
tion rate and the SN rate [1^, [T^I , the prediction for the 
SN neutrino spectra has gone through an evolution over 
the years. Earlier considerations of matter induced res- 
onances were followed by incorporating the 'collective' 
effects due to interaction amongst the neutrinos them- 
selves in the high density central regions of the core. The 
first study of the effect of the collective flavor oscillations 
on DSNB fluxes and the corresponding predicted num- 
ber of events in terrestrial detectors were carried out in 
[g^ where it was demonstrated that the event rate gets 
substantially modified by collective effects. The results 
also showed that observation of the DSNB fluxes at earth 
could shed light on the neutrino mass hierarchy. How- 
ever, since then substantial progress has been made in 
the understanding of the collective effects. In particular, 
it is now clear that the neutrino and antineutrino survival 
probability and hence different split patterns depend cru- 
cially on the relative luminosities of the initial neutrino 
fluxes produced inside the exploding star [36l - [38j . There- 
fore, one can predict the flnal neutrino and antineutrino 
spectra from a given SN with reasonable accuracy only if 
one already has access to the initial flux conditions. This 
complication is further compounded for the DSNB, as 
the DSNB flux comes from a superposition of the fluxes 
from all past SNe. Since the initial flux conditions are 
expected to be sensitive to the properties of the progen- 
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itor star and since we have a whole distribution of stars 
which end up being a SN, it is a comphcated business to 
accurately estimate the DSNB spectra after accounting 
for the collective effects, which are bound to happen in 
almost every SN. 

In this work we incorporate the observation of different 
split patterns in the spectra for the calculation of DSNB 
and do not take the relative (anti)neutrino fluxes to have 
fixed values. The main focus of this work is to check the 
effect of a distribution of supcrnovae with initial flux on 
the measurement of the neutrino mass hierarchy via the 
observation of the DSNB signal. Since the distribution 
of the initial fluxes over all past SNe are not available 
to us, we parametrize this by a log-normal distribution. 
The log-normal distribution has two parameters which 
define the mean and width of the distribution. Since 
they are also unknown, we choose various plausible values 
for them. We calculate the DSNB event rate averaged 
over these distributions. We study how the hierarchy 
measurement is affected when one takes the distribution 
of initial relative fluxes into account and find situations 
where the hierarchy determination may be possible. 



II. THE DIFFUSE SUPERNOVA NEUTRINO 
BACKGROUND 

The differential number flux of DSNB is 



^0 
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where = (l + z)^^_E is the redshiftcd neutrino 
energy observed at earth while E is the neutrino en- 
ergy produced at the source, En is the neutrino flux 
for each core collapse SN, i?sN(^) the cosmic SN rate 
at redshift z, and the Hubble constant taken as i/o = 
70 /i7okm s^i Mpc"^ For the standard A-CDM cosmol- 
ogy, we have matter and dark energy density = 0.27 
and r^A = 0.73, respectively [tH. As Eq. ([T]) suggests 
the DSNB flux at earth depends on two factors: (i) the 
cosmic SN rate and (ii) the initial SN neutrino spectrum 
from each SN. 

The cosmic SN rate is related to the star formation 
rate Rsf{z), through a suitable choice of Initial Mass 
Function (IMF) as Rsn{z)=0.0132x Rsf{z)M-^^,^. 
The IMF takes into account that only stars with masses 
larger than 8Mq result in supernova explosion. For the 
cosmic star formation rate per co-moving volume we take 



i?SF(z) = 0.32 fsNhro- 



„3.4z 



^17^(1 + 2)3 + »A 



,3.82 



+ 45 



(I + z)3/2 



,(2) 



where fg^ is normalization of the order of unity and 
Rsf{z) is in units of Moyr-^Mpc"^ [H,!!!]- The initial 
SN neutrino spectrum emitted from the neutrinosphere 



is parametrized in the form [75| 
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where 0° is the total initial flux estimated for the ini- 
tial luminosity and average energy (£^). The spec- 
tral shape also depends on the energy distribution ip{E), 
which is parametrized by the pinching parameter a. 
In this_ study we use Eu^= 12 McV, i?p^= 15 McV, 
Eu^^ Evy=l8 MeV with Uu^^a^^ = a^^^a^^ = 4 and 
aiy^=Qp^ = 3. Here is a linear combination of Vf^^Vr 
and Vy is the combination orthogonal to Vx] in our case 
Vx and Vy has same flux hence E^^ =E^^ . The average 
energies of the different flux types will also vary from SN 
to SN. However for simplicity, in this work we choose to 
keep the average energies fixed. We assume that 3 x 10^^ 
erg of energy is released in (anti)neutrinos by all SNe. 

The emitted spectrum E^{E) is processed by collec- 
tive flavor oscillation and MSW oscillation effects over 
the huge drop of matter density inside SN. The collective 
oscillations are over within a few 100 km from neutri- 
nosphere whereas the MSW oscillation takes place in the 
region 10^ — 10^ km [tI] for the solar and atmospheric 
mass squared differences. As the collective and MSW 
oscillations are widely separated in space, they can be 
considered independent of each other [35| . Thus the flux 
reaching the MSW resonance region already has the ef- 
fects of the collective oscillations. This assumption may 
not hold in SN models [l^ , where at late times the MSW 
resonance and the collective effects can be simultaneous 
as there the matter density falls substantially compared 
to the early time. In the time integrated DSNB flux it can 
give rise to some corrections. However we have ignored 
these effects as in this work wc followed matter profile 
from studies [s^ with larger matter density [t^ [tJ find- 
ing the two oscillations regimes to be mutually exclusive 
even at late times. 

It has been seen that collective oscillations can give rise 
to different split patterns of the neutrino spectra depend- 
ing on the initial relative flux of i^e and with respect 

to flavor Vx or Vy, so we deflne <Pl,^=^^ and 4'p^ = '^ 
as measures of the relative fluxes [s^. The electron an- 
tineutrino flux beyond the collective region can swap to 
X flavor above some energy (single split) or can swap 
in some energy interval (double split) or even can remain 
unchanged (n o sp lit) depending on the initial relative flux 
(1)1^ and (j>l^ [M^. DSNB is affected differently with 
these different oscillation scenarios. To incorporate the 
effect of collective oscillations we work in a effective two 
flavor scenario with single angle approximation ^. 



^ Multi-angle effects can give rise to kinematical decoherence 
among angular modes and smear the spectral splits [TtII . how- 
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Recent papers [s^ have explored the effect of three 
flavors on the outcome of the split patterns in collective 
oscillations. The three flavor results differ a bit from the 
two flavors only for IH and that too in a small region of 
the initial flux parameter space {(jf^ ), where single 
split (in 3 flavor) appears instead of the double splits (in 
2 flavor). However the observed single split in NH for 
this region remains unchanged in the 3 flavor treatment 
fi^ . Thus in a small part [33, HI] of the parameter space 
((/)J)^,(/)p^), where 2 and 3 flavor results differ, both IH 
and NH have similar split pattern (single split) for the 3 
flavor evolution. In other parts of the initial flux param- 
eter space the split patterns remain the same. Therefore 
while averaging over the whole parameter space the cor- 
rections in the 3 flavor treatment compared to the 2 flavor 
one coming from the small region would not be apprecia- 
ble. The main conclusions of this paper - as we later see 
- therefore remain largely independent of these small cor- 
rections. Moreover the three flavor effects are also found 
to be sensitive to the neutrino-neutrino interaction po- 
tential, the evolution of a three flavor system effectively 
behaves like a two flavor one with the reduction of the 
neutrino-neutrino interaction potential by only one order 
p^ . Hence for such a smaller neutrino- neutrino poten- 
tial the two flavor treatment of the collective evolution is 
reasonable and for larger values the two flavor treatment 
can be again be considered as a reasonable approxima- 
tion when averaging over the whole initial flux parameter 
space is taken into account. 

In this effective two flavor treatment the collective os- 
cillations involve only two flavors {v^, v^), while the other 
flavor [vy) does not evolve. Thus the Flux (F^^ ) after the 
collective oscillations are given by 

= PcFl + (1 - Pc)Fl , Fl = P,Fl + (I - P,)Fl . 

Similarly the quantities F^ and F^ can be estimated 
from the relations 

F'.^ + F-^^ = FO + Fl and Fl + Fl = F^ + F^ . 

As for the other flavor Vy{vy) there is no collective evo- 
lution, hence F,'; ,- ^ = ^ . The quantities Pc and 

Pc are the neutrino and antineutrino survival probability 
after the collective effect, respectively. The only way Vy 
can affect the final neutrino spectrum is by MSW transi- 
tion. As we consider self-induced neutrino oscillations to 
happen independent of the MSW, so the pre-processed 
flux F^^ and the unchanged 'y' flavor will undergo the 
traditional MSW conversions. 

In NH the MSW resonances affect the flux, while 
the Ve flux remains almost unaffected. However for IH, 
MSW affect the flux, and not the flux. Then these 
neutrinos get redshifted while traveling independently as 



ever for spherical symmetry the single angle approximation i,c 
neutrino-neutrino interactions averaged along a single trajectory 
comes out to be a fine approximation as the multi angle deco- 
herence in such a case is weak against the collective features 



Normal hierarchy 

+s12{1-Ph){FS^-FI) + FI. 
F,^ = chPci^:^ , <t>k , E) {Fl - Fl ) + Fl . 



Inverted hierarchy 

F^^=si2P,{(b:^,cl>l^,E){Fl^Fl) + Fl. 

F,, = chPc{4>l,ckk,E)i2PH - l){Fl - Fl) 
+cU1-Ph){FI-FI) + FI. 

TABLE I: Electron neutrino and antineutrino spectra emerg- 
ing from a SN. 



mass-eigenstates until they reach earth, here they are 
detected as flavor eigenstates before or after having un- 
dergone regeneration inside the earth. The fluxes {F,^^ 
and Fp_.) arriving at earth after both the collective and 
MSW oscillation are given in Table I. 

In Table I, Ph is the effective jump probability be- 
tween the neutrino mass eigenstates due to the atmo- 
spheric mass squared driven MSW resonance. It takes 
a value between and 1 depending on the value of the 
mixing angle ^13. For ^13 large (sin^ ^13 ^ 0.01) Ph — 0, 
while for smafl 6I13 (sin^ ^13 ^ 10"^) Pj? ~ 1 (zl. For 
the smaller limit of ^13 the effective jump probability 
has a non-trivial dependence on energy and time, due 
to multiple resonances [14, 19] and turbulence [16, 17, 
18]. Since these effects occur in the small time- window 
when the shock wave is in the resonance region one can 
neglect this sub-leading effect in the time integrated ob- 
servablcs like DSNB events. The quantities sfj and cfj 
stand for sin^ 6*12 (taken to be 0.3 for numerical studies) 
and cos^ ^12, respectively. The neutrino and antineutrino 
survival probability after the collective effect are calcu- 
lated numerically taking Am^ = 3 x 10"'^ eV^ and a 
small effective mixing angle of 10^^. However it should 
be noted that the collective probabilities Pc and Pc are 
almost independent of the specific non zero values of the 
mixing angle [l^, H^, [s^ . Thus the effective mixing an- 
gle in the matter for both the limits sin^ 613 ^ 0.01 and 
sin^ 613 ^ 10^® will give rise to similar collective fea- 
tures. We checked explicitly that the results remain the 
same when one varies the effective mixing angle by or- 
ders of magnitude. Hence for the collective evolutions 
we consider a representative value of 10~^ as the effec- 
tive mixing angle for both the limits. For further details 
of analysis of the collective effects we refer the reader to 
[sst . ^ As discussed before, Pc and/or Pc, which is a 
function of the neutrino energy E, show pattern of sud- 
den change between and 1 , leading to sudden change in 
the neutrino and/or antineutrino spectra. Depending on 



^ Note that Pc = Pc = 1 gives back the SN flux without collective 
oscillation Q| 
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FIG. 1: The four specimen distributions used in the calculations. 
The y-axis gives the distribution of supernova in arbitrary units. 



the number of times the value of Pc (and/or Pc) changes, 
we can have more than one sudden swap of the neutrino 
(and/or antineutrino) spectra, referred to in the htera- 
ture as multiple split |36l438j . The split patterns are now 
known to be crucially dependent on the initial relative 
flux densities of the neutrino and antineutrino. The ini- 
tial relative neutrino and antineutrino flux densities are 
expected to vary among different SN. 

This variation in (/)^ and 0p in different SNe might 
have origin in difference of progenitor mass, different lu- 
minosity etc. In fact even different simulations allow wide 
variation of (/)^ and (j)^ within the 2 fold uncertainty 

around equipartition ( < (ff^,^ < < </>p^ < 

•) [ll- Most models predict l^,^ ~ /p^, where 

and Ij}^ are the relative luminosity of i^e and Df, 

respectively (STj . However, the combined luminosity of 
z/^, Vfj^ and is seen to be rather disparate between 
the different model results'^. The most reliable way to re- 
construct the relative luminosity distribution function in 
principle would be from direct observation of SN events 
along with their neutrino signal. However, as yet only 
SN1987A has been observed along with the detection of 
its neutrinos/antineutrinos. We require to know neutrino 
fluxes of different flavors from a number of galactic SN 
with a range of stellar mass and initial conditions before 
collapse to have information about the possible variation 
in (pl, and 0p . This might take decades and is clearly 
not possible in the near future. So we propose in this 
paper that the SN events contributing to DSNB have 
various different values of 0^ and 0p . For quantitative 
estimates we assume specific distributions for them. We 



^ See for e.g. the compilation of model results in [80| . 



take (^J^^ and i^p^ distributed log-normally, defined by the 
parameters (/xi,(Ti) and {^12,(^2) i-e, 

D^{^i2,(^2) = jK^—r • (4) 

\/2t:(J2<Pu^ 

We choose a range of values for ji such that the ex- 
pectation (mean) values of 0^ and (jf^^ are compatible 
with either Lawrence Livermore (equipartition) or Garch- 
ing simulations. The parameter a determines the width 
of the distribution. Since the variation of (pl, and 0p 
is expected to be within the 2 fold uncertainty around 
equipartition Q, the cr's for the distributions are chosen 
such that the distribution is well within the 2 fold un- 
certainty of the expectation value. In order to show the 
effect of the choice of the distribution function we simu- 
late our results for four widely different distributions for 
(ff^, and . They are chosen as follows: 

1. The expectation is the equipartition value, i.e 
1.5 for (jfjj and 1.2 for i^p and the variation is 
within 2 fold uncertainty of 1.5 and 1.2 respec- 
tively. The corresponding values of (/i, u) turn out 
to be (0.39,0.21) for neutrinos and (0.16,0.24) for 
antineutrinos. We denote this distributions by D], 
and D\^. 

2. The expectation is same for both (jf^ and (/>p and is 
taken as 1, the a is chosen as stated above. Hence 
the distributions for v^. and v^. are identical and is 
defined by the parameters (/.t, cr) = (—0.03,0.28). 
We denote this distribution by and Z?? . 

3. The expectation is taken to be the same as the 
Garching simulations [zHl, i,e 0.8 for both the 
i/g and and the distributions parametrized by 
(/x,ct) = (—0.23,0.20) are the same. We denote 
this distribution by D^^ and D^^ . 

4. We consider the distribution to be constant be- 
tween 0.45-2.55 for and 0.45-2.05 for pl^, with 
a normalization factor of and j-q respectively. 
We denote this case by and . 

The four specimen distributions for and are 
shown in Fig. [TJ By definition the log-normal distribu- 
tion functions , and arc normalized to unity. 
The choice of our normalization factors for the uniform 
distribution functions in D"* ensures that they are nor- 
malized to unity as well. We can see from the figure that 
D'^ has the narrowest spread in the relative fiux while the 
uniform distribution has the widest. For simplicity we 
have chosen distribution functions that are independent 
of the redshift z. 
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The differential number flux of DSNB with the i*^ dis- 
tributions is given by 



c 



,(5) 



where Dl, {(f>l, ) and Dp (0p ) is the number of super- 
novae in between the interval 0[) to t/)^ -I- and 0p 
to , respectively. {E, (jf^^ , (jf^^ ) is the neutrino 

flux of each SN with initial relative flux combination 0^ 
and at a redshift z. Zmax is taken as 5, where as 
4>amm ^^'^ 4>amax S'^'^ chosen from the 2 fold uncertainty 
around expectation of the respective distribution. The 
upper panel in Fig. [5] shows the antineutrino DSNB flux 
(in logarithmic scale) for NH and the lower panel shows 
the flux for IH with Ph= 1. In both panels flux is shown 
for all of our different distribution (D^-D^) of (p^ and 
. In Fig. [5] we also plotted flux for the case without 
any distribution of and 0p but speciflc value (0.8, 
0.8) for the initial relative fluxes (6^ ,6r, )■ From the 
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FIG. 2: Ve fluxes for NH(upper panel) and IH (lower panel) 
with small dis, i.e Ph= 1. 



Model 


Hierarchy 


2.5 Mton-yr 
(19.3 - 30.0) 
(MeV) 


GD-f 2.5 Mton-yr 
(10.0 - 30.0) 
(MeV) 




NH 


76±9 


280±17 


IH {Ph = 0) 


92±10 


319±18 


IH {Ph = 1) 


80±9 


288±17 


(1.5,1.2) 


NH 


75±9 


281±17 


IH {Ph = 0) 


98±10 


332±19 


IH {Ph = 1) 


75±9 


280±17 




NH 


88±10 


307±18 


IH {Ph = 0) 


103±11 


350±19 


IH {Ph = 1) 


76±9 


280±17 


(1.0,1.0) 


NH 


102±11 


340±19 


IH {Ph = 0) 


103±11 


341±19 


IH {Ph = 1) 


77±9 


295±18 




NH 


98±10 


334±19 


IH {Ph = 0) 


113±11 


379±20 


IH {Ph = 1) 


70±9 


258±16 


(0.8,0.8) 


NH 


112±11 


378±20 


IH {Ph = 0) 


113±11 


378±20 


IH {Ph = 1) 


70±9 


260±17 




NH 


82±9 


297±18 


IH {Ph = 0) 


97±10 


335±19 


IH {Ph = 1) 


79±9 


283±17 


(1.5,1.25) 


NH 


77±9 


286±18 


IH {Ph = 0) 


97±10 


330±19 


IH {Ph = 1) 


76±9 


285±17 



TABLE II: Number of expected events per 2.5 Megaton-year 
of a water Cherenkov with SK like resolution and in a similar 
detector with Gadolinium loaded. 



flux figures it is evident that the different distributions 
give similar DSNB flux. To check whether the small dif- 
ferences in their profile are measurable or not we next 
calculate the corresponding event rate in different detec- 
tors. 



III. DSNB EVENT RATES AND HIERARCHY 
MEASUREMENT 

Since the DSNB fluxes are very small, they are ex- 
pected to be observed in either very large detectors or in 
reasonable size detectors with very large exposure times. 
The only reasonable size detector running currently is 
the Super-Kamiokande (SK) [g^]. Amongst the proposed 
large detectors which have low ener gy t hreshold are the 

S3, Gadolinium 



81 



8J] (which could be 



megaton water Cherenkov detectors 
enriched water Cherenkov detectors 
SK or possibly an even larger water detector), very large 
liquid scintillator detector (LENA) [H^ or very large liq- 
uid Argon detector [ssl - ISTj . Apart from the liquid Argon 
detector, none of the other detectors hold much promise 
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2.25 



FIG. 3: Upper panel shows the number of expected DSNB events 
in 2.5 MtonYr of Gadohnium loaded water detector for normal 
(NH) and inverted (IH) hierarchies (for IH, Pfj = 1 ) as a function 
°f 'Pit {= )• The lower panel shows the change of r; as a function 

of «J= 



for the detection of DSNB Ve- Therefore, in rest of this 
discussion we will focus only on the detection possibilities 
for the DSNB fluxes. All detector technologies men- 
tioned above use the inverse beta decay v^+p — > n-\-e^ 
interaction for detecting DSNB fluxes. The only differ- 
ence between them would be in terms of their energy 
window of sensitivity for the DSNB. Water detectors use 
the neutrino energy {E^) window 19.3 MeV to 30.0 MeV 
[67j . whereas for Gadolinium loaded water detectors the 
detection window is between 10 MeV and 30 MeV js^ . 
Liquid scintillator on the other hand has a sensitivity 
range between 10 MeV and 25 MeV We show the 
total number of events for 2.5 Mton-yr exposure in wa- 
ter detectors and 2.5 Mton-yr exposure in Gadolinium 
loaded water detectors in the third and fourth columns 
of Table II. The number of events are shown for the four 
different specimen distributions. Results are shown sep- 
arately for Ph ~ 1 and Ph = when the neutrino mass 
hierarchy is inverted (IH) . For the NH, there is no de- 
pendence of the Pe flux on Ph, as one can see from Table 
I. In each case we also show for comparison the number 
of events expected when (p^^ and (pp are fixed at their 
respective mean values of the distribution ^. Observing 
a variation of the number of events changing the and 
for these examples without any distribution is an 
important indicator of how inclusion of distribution can 
change the DSNB event numbers. 

We next turn our attention to the possibility of mea- 
suring the neutrino mass hierarchy using DSNB detec- 



tion. For the distribution examples D^, and Z?* the 
difference of number of events between NH and IH (both 
Ph = 1 and Ph ~ 0) seems to be within statistical un- 
certainty. Whereas for the distribution D^, though the 
difference in number of events between NH and IH (with 
Ph = 1) decreases significantly compared to the without 
distribution case, the difference is still greater than the 
statistical uncertainty. To probe this variation of hierar- 
chy sensitivity with the distribution further let us define 
the quantity 



NH 



Nih\ 



Nnh 



(6) 



which gives a measure of the hierarchy sensitivity of the 
experiment. The quantities Nnh and Nih are the num- 
ber of expected DSNB events when the hierarchy is nor- 
mal and inverted, respectively. The quantity rj is defi- 
nitely not an observable as it depends on both N^h and 
NjH ■ However it helps describing the variation of hierar- 
chy sensitivity of DSNB in a systematic manner. There 
are two ingredients of interest in Table II. Firstly, we can 
see that the relative hierarchy difference rj depends on 
the mean value of (f>l, and 0^ . Secondly, for a given 
< 4>l, > and < > it also depends on the distribution 
function involved. 

For a better understanding of the first issue, we show 
in the upper panel of Fig. [3] the variation of the number 
of events for both the NH and IH cases, as a function of 
the relative luminosity factor. The analysis is again done 
for the small mixing angle limit {Ph = 1) for IH, as it is 
the more challenging limit from the experimental point 
of view. We show this for 2.5 McgatonYcars of Gadolin- 
ium doped water detector^. For simplicity we have taken 
= in this figure and we do not take any distri- 
bution function into account. The lower panel shows the 
corresponding 77 as a function of c/)^^ = . We see that 77 
has a very complicated dependence on the relative initial 
luminosity functions. The variation of 77 with 0^s ac- 
tually depends on several factors like difference between 
split pattern of NH and IH, the split/swap energies, the 
initial relative flux of u^, 9^ and Vx- It is rather high for 
very low values of and (j)^^. It starts decreasing as 
the value of (f)^, and (j)^ increase until it becomes zero 
around (j)"^^ = ~ 1.05, thereafter it increases for a 
short while until it reaches a (local) maximum at around 
4>l, = 4>^ — 1.25. Beyond that the value of rj decreases 
again reaching a second zero at around = ~ 1.5. 
After that it increases monotonically. 

The most noteworthy thing in the upper panel of this 
figure is that the number of events for inverted hierarchy 



* For the uniform distribution we took the central values of the 
widths as their representative mean. 



^ Most of the results shown in this paper is for megaton class 
Gadolinium doped water detector to show the impact of the 
relative initial luminosity and its distribution on the hierarchy 
measurement using DSNB fluxes. The corresponding sensitivi- 
ties for smaller scale detectors can be calculated trivially using 
these numbers. 
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FIG. 4: Relative difference {r]) of number of expected events in 
normal and inverted hierarchy, per 2.5 Megaton-year in a Gadolin- 
ium loaded SK like detector. For the IH case we consider small 
mixing angle i,e Ph = 1. The above panels show r] for four possi- 
ble cases with a being 0.00, 0.05, 0.10, 0.15 and 0.20. The x axis 
denotes the expectation of 0^ and for the chosen sigma. Here 
expectation of ip^ and (f>^ are taken to be same i,e fii = fi2. The 
errors shown are the statistical errors only. 



FIG. 5: Plot of 



\Nkh-Nih\ 



with < 



>(= < 't'l, » for 



o'=0.20. The ratio is shown for a Gadolinium loaded SK like de- 
tector with different values of the exposure i.e. 2.5, 5.0, 7.5 and 
10.0 Mton-Yr. Horizontal line at ratio 1 denotes that statistical 
error is of the same order as the difference between two hierarchies. 
Higher values like 2, 3 will be of more prominence as other uncer- 
tainties from flux, luminosities and detector systematics should be 
considered alongwith statistical error . 



has almost linear dependence on (j)^^ = 0^ on both sides 
of the maximum which comes around = '/'p — 1.25. 
For the normal hierarchy one sees departure from linear- 
ity around = 0p ~ 1.1. This feature is crucial in 
determining the effect of the distribution function on the 
hierarchy sensitivity. The effect of taking the distribution 
function into account boils down to creating a weighted 
average of the number of events, where the weights are 
determined by the distribution itself. For the log-normal 
case the weights are driven by the mean and width of 
the distribution, which we parametrize in terms of fi and 
a. The effect of any distribution can thus be understood 
with the help of Fig. [3] 

To show the effect of the distribution function we con- 
tinue to stick to the simplified scenario where = 0p 
and show in Fig. |4]the relative difference as a function 
of the expectation value (0^ )(= (0^ )). Here again we 
take Ph = 1 for IH. We consider log-normal distribution 
for the relative luminosities ^ and show 7] for five different 
values of a which controls the width of the distribution. 
The values and error bars on rj correspond to a statistics 
of 2.5 Megaton- Year data in Gadolinium loaded water 
detector. The lowest panel with a = corresponds to 
the case where we keep (0^^ ) fixed and for this case 
there is no effect of the distribution. This case is similar 



Our conclusions remain fairly robust against the choice of the 
distribution function. We have explicitly checked this by repeat- 
ing Fig. |4] with uniform distribution and normal distribution. 
However, we do not present those results as they follow the same 
pattern that we get for the log-normal distribution. 



to that in the lower panel of Fig. |3l We see that without 
the effect of the distribution function almost all values 
of (j)^^ (= ) would give hierarchy sensitivity to at least 
la C.L., while for lower values of the relative luminosity 
the sensitivity can be seen to be rather good. As we in- 
crease (J the sensitivity is seen to go down for all value 
of the relative luminosity. We find that even for very 
small values of tr = 0.05, there is almost no hierarchy 
sensitivity for (0^ )(= (0^ )) ^1. As cr increases this 
pattern remains the same, though the sensitivity keeps 
falling for all values of ) (= (0p )). In the above anal- 
ysis, a is considered in a small range. The idea was to 
avoid deviating too much from the simulated values of 
the relative initial luminosities. If the actual variations 
of 0^s for all past supernovae are much larger than the 
range considered here, then the difference between pre- 
dicted events for NH and IH would decrease even further, 
and this could washout the hierarchy sensitivity even for 
the low (pl, cases. 

Finally in Fig. [5] we display how we can make sure 
quantitatively that the observed number of events be- 
long to one of the hierarchies and not the other one and 
see its dependence on the exposure in Mton-Yr. Fig. [5] 
plots the ratio for four different values of 

the exposure i.e. 2.5, 5.0, 7.5 and 10.0 Mton-Yr as a 
function of the expectation values of the relative flux, 
< (f>l^ > taken as equal to < (j)^^ >. This is for the 
case Ph=1 and with (7=0.20. This shows that with 2.5 
Mton-Yr, the relative flux expectations with values be- 
low 0.95 (ratio ~ 1) can distinguish NH from IH (at least 
in principle) as the statistical error is smaller than the 
number of event difference. However to achieve a better 
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sensitivity i.e ratio 2 the upper limit for the relative 
flux expectations is about 0.77. For ratio ^ 3 the value 
is even lower. However with larger values of exposures 
like 5.0 or 7.5 Mton-Yr one can still make a distinction 
between hierarchies when the centroid is below 0.75 and 
0.82 respectively. But it is clear from the figure that by 
increasing the Mton-Yr one gets only a slow increase in 
the allowed parameter space. 

IV. CONCLUSION 

In this work we studied the prospects of measuring the 
neutrino mass hierarchy from observation of the DSNB 
signal in terrestrial detectors. While such studies have 
been performed before by different groups including ours, 
our work is unique as this is the first time that distribu- 
tion of the source SN with initial relative neutrino and 
antincutrino fluxes has been taken into account. It is 
natural that different SN would emit neutrino and an- 
tincutrino fluxes with slightly different initial conditions, 
depending on the properties of the progenitor star. This 
is particularly relevant in the context of collective oscilla- 
tions, where the multiple split patterns depend crucially 
on the initial relative fluxes. Since the actual distribu- 
tion function of SN with the initial relative fluxes are 
unknown, we chose four specimen distribution functions, 
which have a mean corresponding to the value from SN 
simulations and a width such that almost all the val- 
ues are within a factor of two of the mean value. We 
worked with three log-normal and one uniform distri- 
bution. We presented the DSNB fluxes for all the four 
distributions for both normal and inverted hierarchies. 
We calculated the total predicted number of De events 
in water detectors, both with and without Gadolinium. 
The log-normal distribution is characterized by its mean 
value and its variance. These are parametrized in terms 
of the variables n and a. We studied the dependence 
of the hierarchy sensitivity to the mean and variance of 
the log-normal distribution function. We concluded that 
the hierarchy sensitivity in this experiment had a cru- 
cial dependence on the mean value of the relative initial 
luminosity 0^ and 0^ . The sensitivity has a predomi- 
nantly non-linear dependence on ) and (0^ ), being 
higher for lower values of these quantities. The effect of 
the variance parametrized by tr is to reduce the hierarchy 
sensitivity for all values of the mean {cj)1J and {(jy^J. We 



found that even for very moderate values of cr ~ 0.05, 
there is almost no hierarchy sensitivity in the very small 
mixing angle limit for {4>U = ((/>pj ^ 1. 

Finally we discuss the issue of using an effective two 
flavor evolution followed in this paper. As mentioned 
earlier if the analysis is done using all three flavors only 
the 'double splits' for IH arising in a small part of the 
((^^ , (j)^ ) parameter space seen in the two flavor case 
change to 'single split' while things remain unchanged 
for NH. This effect would have been important for calcu- 
lations for a single SN event if the initial flux ratios fell 
in this part of the parameter space. However for DSNB 
with averaging done over the SN distributions (Di) the 
three flavor treatment will result in small corrections to 
the effective two flavor results presented here. 

To summarize, using some reasonable simplifying as- 
sumptions we for the first time take into account, for 
calculating DSNB, a distribution in the relative fluxes in 
SN neutrinos and study its effect on the ability to distin- 
guish the two hierarchies through a future observation 
of DSNB. A more detailed and rigorous study involv- 
ing multi-angle effects and three flavors will be able to 
give more exhaustive answers to many of the aspects and 
should be persued in future. But the model calculation 
carried out in this paper asserts that for DSNB which is 
a blend of many SN, not all identical, it is more difficult 
to distinguish IH from NH compared to a single SN or 
the hypothetical case of DSNB made of supernovae with 
identical flux ratios. A few favorable scenarios for this 
hierarchy distinguishability are indicated. We feel that 
the central finding of this first study is robust and will 
survive more detailed analysis 
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